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A commonly used theoretical definition of superfluidity in the ground state of a Bose gas is based 
on the response of the system to an imposed velocity field or, equivalently, to twisted boundary 
conditions in a box. We are able to carry out this program in the case of a dilute interacting Bose 
gas in a trap, and we prove that a gas with repulsive interactions is f00% superfluid in the dilute 
limit in which the Gross-Pitaevskii equation is exact. This is the first example in an experimentally 
realistic continuum model in which superfluidity is rigorously verified. 



I. INTRODUCTION 

The phenomenological two-fluid model of superfluidity 
(see, e.g.Ji) is based on the idea that the particle density 
p is composed of two parts, the density p s of the inviscid 
superfluid and the normal fluid density p n . If an exter- 
nal velocity field is imposed on the fluid (for instance by 
moving the walls of the container) only the viscous nor- 
mal component responds to the velocity field, while the 
superfluid component stays at rest. In accord with these 
ideas the superfluid density in the ground state is often 
defined as followscl: Let Eq denote the ground state en- 
ergy of the system in the rest frame and E' Q the ground 
state energy, measured in the moving frame, when a ve- 
locity field v is imposed. Then for small v 



f = ^ + Wp)^v 2 + 0(|v| 



(1) 



where N is the particle number and m the particle mass. 
At positive temperatures the ground state energy should 
be replaced by the free energy. (Remark: It is important 
here that ([j]) holds uniformly for all large N; i.e., that the 
error term 0(|v| ) can be bounded independently of N. 
For fixed N and a finite box, Eq. ([!]) with p s /p = 1 al- 
ways holds for a Bose gas with an arbitrary interaction if 
v is small-enough, owing to the discreteness of the energy 
spectrumd.) There are other definitions of the superfluid 
density that may lead to slightly different resultso, but 
this is the one we shall use in this paper. We shall not 
dwell on this issue here, since it is not clear that there 
is a "one-size-fits-all" definition of superfluidity. For in- 
stance, in the definition we use here the ideal Bose gas 
is a perfect superfluid in its ground state, whereas the 
definition of Landau in terms of a linear dispersion rela- 
tion of elementary excitations would indicate otherwise. 
We emphasize that we are not advocating any particular 
approach to the superfluidity question; our contribution 
here consists in taking one standard definition and mak- 
ing its consequences explicit. 

One of the unresolved issues in the theory of superfluid- 
ity is its relation to Bose-Einstein condensation (BEC). 
It has been argued that in general neither condition is 



necessary for the other (c.f., e.g.|ffl'Ll). A simple exam- 
ple illustrating the fact that BEC is not necessary for 
superfluidity is the ID hard-core Bose gas. This system 
is well known to have a spectrum like that of an ideal 
Fermi gasu, and it is easy to see that it is superfluid in 
its ground staJje-in the sense of ([j]). On the other hand, 
it has no BECoo. The definition of the superfluid veloc- 
ity as th£j_gradient of the phase of the condensate wave 
functionBEil is clearly not applicable in such cases. 

We do not give a historical overview of superfluid- 
ity because excellent review articles are availableOE3. 
While the early investigations of superfluidity and Bose- 
Einstein condensation were mostly concerned with liq- 
uid Helium 4, it has become possible in recent years 
to study these phenomena in dilute trapped gases of 
alkali atomscJ. The experimental success in realizing 
BEC in such gases has led to a large number of theo- 
retical papers on this subject. Most of these works take 
BEC for granted and start off with the Gross-Pitaevskii 
(GP) equation to describe the condensate wave function. 
A rigorous justification of these assumptions is however 
a difficult task, and only very recently BEC has been 
rigorously proved for a physically realistic many-body 
Hamiltoniaro. It is clearly of interest to show that su- 
perfluidity also holds in this model and this is what we ac- 
complish here. We prove that the ground state of a Bose 
gas with short range, repulsive interaction is 100% su- 
perfluid in the dilute limit in which the Gross-Pitaevskii 
description of the gas is exact. This is the limit in which 
the particle number tends to infinity, but the ratio Na/L, 
where a is the scattering length of the interaction poten- 
tial and L the box size, is kept fixed. (The significance of 
the parameter N a /L is that it is the ratio of the ground 
state energy per particle, ~ Na/L 3 , to the lowest exci- 
tation energy in the box, ~ 1/L 2 .) In addition we show 
that the gas remains 100% Bose-Einstein condensed in 
this limit, also for a finite velocity v. Both results can 
be generalized from periodic boxes to (non-constant) ve- 
locity fields in a cylindrical geometry. 

The results of this paper have been conjectured for 
many years, and it is gratifying that they can be proved 
from first principles. They represent the first example of 
a rigorous verification of superfluidity in an experimen- 
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tally realistic continuum model. 

We wish to emphasize that in this GP limit the fact 
that there is 100% condensation does not mean that no 
significant interactions occur. The kinetic and poten- 
tial energies can differ markedly from that obtained with 
a simple variational function that is an A-fold product 
of one-body condensate wave functions. This assertion 
might seem paradoxical, and the explanation is that near 
the GP limit the region in which the wave function differs 
from the condensate function has a tiny volume that goes 
to zero as N — > oo. Nevertheless, the interaction energy, 
which is proportional to N, resides in this tiny volume. 



L. This means that a should vary like 1/N which can 



be achieved by writing v(r) 



~Vi(r/a), where «i is a 



fixed potential with scattering length 1, while a changes 
with N. 

Theorem 1 (Superfluidity of homogeneous gas). 

For \ip\ < 7r 



lim 



Eo(N,a,<p) 
N 



An flap + ji 



L 2 



(6) 



in the limit N — > oo with Na/L and L fixed. Here p — 
N/L 3 , so ap is fixed too. In the same limit, for \cp\ < tt, 



II. SETTING AND MAIN RESULTS 



We consider a Bose gas with the Hamiltonian 



N 



H 



N 



E 

l<i<j<JV 



(2) 



where /i = h 2 /(2m) and the interaction potential v is 
nonnegative and of finite range. The two-body scat- 
tering length of v is denoted by a. The Hamiltonian 
acts on totally symmetric functions W of N variables 
Yi = (xi,yi,Zi) e JC C R 3 , where /C denotes the cube 
[0, L] 3 of side length L. (We could easily use a cuboid of 
sides Li,L2,i3 instead.) We assume periodic boundary 
conditions in all three coordinate directions. 

Imposing an external velocity field v = (0,0, ±|v|) 
means that the momentum operator p = — ifi,V is re- 
placed by by p — mv, retaining the periodic boundary 
conditions. The Hamiltonian in the moving frame is thus 



N 



-.^(Vj+itp/L) + E v (\^~ r J 



l<i<j<N 



(3) 



where tp = (0, 0, ip) and the dimensionless phase p is 
connected to the velocity v by 



±|v|Lm 



(4) 



Let Eq(N,o, tp) denote the ground state energy of (j|) 
with periodic boundary conditions. Obviously it is no 
restriction to consider only the case — n < tp < ir, since 
Eq is periodic in ip with period 2ir. For the ground 
state of H' N , let 7at be its one-particle reduced density 
matrix 



7jv(r,r') =N / V (r,r 2 , . . . ,r N ) 



x *g(r',r 2 , . . . ,r N )dr 2 ■ ■ -dr N . (5) 

We are interested in the Gross- Pitaevskii (GP) limit 
N — > oo with Na/L fixed. We also fix the box size 



lim -L 7jv ( r y) = ^ (7) 

in trace class norm, i.e., limyv^oo Trace [ ryjv/^V — 
\L- 3 / 2 ){L- 3 ' 2 \ |] =0. 

Note that, by the definition ([!]) of p s and Eq. (^), Eq. 
(||) means that p s — p, i.e., there is,4|00% superfluidity. 
For ip = 0, Eq. (||) was_first proved inta. Eq. (0) for ip = 
is the BEC proved intd. 

Remarks. 1. By a unitary gauge transformation, 



([/*)(n,...,r w ) 



•)/ L *( ri; ..., rjv) , (g) 



the passage from (^) to (^) is equivalent to replacing pe- 
riodic boundary conditions in a box by the twisted bound- 
ary condition 

(ri + (0, 0, L), r 2 , . . . , r N ) = e^*(n, r 2 , ...,r N ) (9) 

in the direction of the velocity field, while retaining the 
original Hamiltionian (Q). 

2. The criterion \ip\ < 7r means that |v| < ■nh/{mL). 
The corresponding energy (jnL)) 2 is the gap in 
the excitation spectrum of the one-particle Hamiltonian 
in the finite-size system. 

3. The reason that we have to restrict ourselves to 
| ip\ < 7r in the second part of Theorem |l| is that for | ip\ = tt 
there are two ground states of the operator (V + itp/L) 2 
with periodic boundary conditions. All we can say in this 
case is that there is a subsequence of jn that converges 
to a density matrix of rank < 2, whose range is spanned 
by these two functions. 

Theorem ^ can be generalized in various ways to a 
physically more realistic setting. As an example, let C 
be a finite cylinder based on an annulus centered at the 
origin. Given a bounded, real function a{r, z) let A be the 
vector field (in polar coordinates) A{r, 9, z) = tpa(r, z)eg, 
where eg is the unit vector in the 9 direction. We also 
allow for a bounded external potential V(r, z) that does 
not depend on 9. _ 

Using the methods of Appendix A irJla, it is not diffi- 
cult to see that there exists a <po > 0, depending only on 
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C and a(r, z), such that for all | p>\ < ipQ there is a unique 
minimizer <j> of the Gross-Pitaevskii functional 

£ GP W= I (H(V + L4(r))<Hr)| 2 

+ y(r)|</>(r)| 2 +47r^7Va|0(r)| 4 )d 3 r (10) 

under the normalization condition j \<fi\ 2 — 1. This min- 
imizer does not depend on 9, and can be chosen to be 
positive, for the following reason: The relevant term in 
the kinetic energy is T = —r~ 2 [d/d9 + itpr a(r, z)] 2 . If 
\ip r a(r, z)\ < 1/2, it is easy to see that T > tp 2 a(r, z) 2 , in 
which case, without raising the energy, we can replace <f> 
by the square root of the-$-average of \<fi\ 2 . This can only 
lower the kinetic energyO and, by convexity of x — » x 2 , 
this also lowers the 4> 4 term. 

We denote the ground state energy of £ GP by E GP , de- 
pending on Na and tp. The following Theorem ^| concerns 
the ground state energy E of 

N 
3=1 

+ «(i r <- r ii). (ii) 

l<i<j<N 

with Neumann boundary conditions on C, and the 
one-particle reduced density matrix of the ground 
state, respectively. Different boundary conditions can 
be treated in the same manner, if they are also used in 
©• 

Remark. As a special case, consider a uniformly rotat- 
ing system. In this case A(r) = ipreg, where 2ip is the 
angular velocity. is the Hamiltonian in the rotating 
frame, but with external potential V{r) + [iA(r) 2 (see 
e.g.y (p. 131)). 

Theorem 2 (Superfluidity in a cylinder). For 

\<P\ < 

lim ^^=E^ { NaM (12) 

N— too 1\ 

in the limit N — * oo with Na fixed. In the same limit, 

lim l 7A ,(r,r')=0 GP (r)^ GP (r') (13) 

N— »oo iV 

in trace class norm, i.e.. limjv^oo Trace [ \Jn/N — 
|0 GP )(0 GP ||] -0. 

In the case of a uniformly rotating system, where 2tp 
is the angular velocity, the condition \ip\ < ipo in partic- 
ular means that the angular velocity-is smaller than the 
critical velocity for creating vorticesEJ. 

Remark. In the special case of the curl-free vector po- 
tential A(r,9) = (pr^ 1 ^, i.e., a(r,z) — r" 1 , one can say 
more about the role of po . In this case, there is a unique 
GP minimizer for all <p ^ Z + | , whereas there are two 
minimizers for ip E Z + h . Part two of Theorem || holds 
in this special case for all ip Z + |, and (|lj) is true 
even for all tp. 



III. PROOFS 

In the following, we will present only a proof of Theo- 
rem |l| for simplicity. Theorem [2] can be proved .using the 
same methods, and additionally the methods oftj to deal 
with the inhomogeneity of the system. 

Before giving the formal proofs, we outlinp_thc main 
ideas. The strategy is related to the one irO, but re- 
quires substantial generalizations of the techniques. A 
crucial element of the proof, stated in Lemma [l] below, 
is the fact that the interaction energy can be localized 
in small balls around each particle. This pact uses a 
Lemma of DysonEj, and its generalization inlij, which 
converts a strong short range potential into a soft poten- 
tial. This Lemma can be also be applied to the case of an 
external velocity field, i.e., a U(l) gauge field in the ki- 
netic term of-the Hamiltonian, owing to the "diamagnetic 
inequaltity"tL3. This inequaltity says that the additional 
gauge field increases the kinetic energy density. 

The second main part of the proof is the generalized 
Poincare inequality given in Lemma ^|. We recall that 
an essential iaeredient of the proof of Bose-Einstein con- 
densation inEJ was showing that the fact that the ki- 
netic energy density is small in most of the configuration 
space implies that the one-body reduced density matrix 
is essentially constant. The difficulty comes from the 
fact that the region in which the kinetic energy is small 
can, in principle, be broken up into disjoint subregions, 
thereby permitting different constants in different subre- 
gions. The fact that this does not happen is the content 
of the generalized Poincare inequality. In the present case 
we have an additional complication coming from the im- 
posed gauge field. The old Poincare inequality does not 
suffice; one now has to measure the kinetic energy den- 
sity relative to the lowest energy of a free particle in the 
gauge field rather than to zero. This is an essential com- 
plication. While the previous (generalized) Poincare in- 
equality could, after some argumentation, be related to 
the standard Poincare inequalityHI, this new one, with 
the gauge field, requires a totally new proof. 

Proof of Theorem As inE! we define Y — (An/ 3) pa 3 . 
Note that in the limit considered, Y ~ ./V~ 2 . We first 
consider the upper bound to Eq. Using the ground state 
for <P — as a trial function, we immediately get 

9 

Eo(N,a,<p)<(* a ,H',f9 )=Eo(N,a,0)+Nii^ , 

(14) 

since (* , Vj* ) = 0. FromE! we know that 
E o (N,a,0) < 4TT^Npa(l + const. Y 1 / 3 ), which has the 
right form as N — * oo. 

For the lower bound to the ground state energy we 
need the following Lemma. 

Lemma 1 (Localization of energy). For all sym- 
metric, normalized wave functions \&(ri, . . . , rjv) with 
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periodic boundary conditions onK, and for N > Y 1 / 17 J 

— (^>,H' N ^) > (1 -const. F 1 / 17 ) (4 W a + 
At/ dX I dri|(Vi+ty>/L)tf(ri,X)| 

where X = (r2, . . . , r^), dX = 11^=2 ^ r j> an d 



(15) 



fix = s ri : min Iri — r,-| > i? 
1 i>2 



(16) 



with R = aY- 5 / 17 . 



Proof. Since '5 is symmetric, the left side of (JTJ]) can be 
written as 



/ dX f dr 1 L\(V 1 +itp/L)*(r 1 ,X) 



i^^-r.-Dl^r^X)! 1 

J">2 



For any e > and R > this is 



> eT + (1 - e)(T in +/) + (!- e)T° 



with 



T = fi[ dX [ dn|Vi|*(ri,X)|| a , 



(17) 



(18) 



(19) 



T in = [i dX dn|Vi|*(n,x)|p , (20) 



The following LemmaJ is needed for a lower bound 
on the second term in Jl5). It is stated for /C the L x 
L x L-cube with periodic boundary conditions, but it 
can be generalized to arbitrary connected sets K, that are 
sufficiently nice so that the Rellich-Kondrashov Theorem 
(seeEJJ (Thm. 8.9)) holds on K. In .particular, this is the 
case if K. has the 'cone property'EZl. Another possible 
generalization is— tp include general bounded vector fields 
replacing tp, seeEil. 

If fi is any subset of JC we shall denote J Q f*(r)g(r)dr 

by (f,9)n and (/,/)^ /2 by \\f\\ L ^(n). We also denote 
V + itp by V v for short. 

Lemma 2 (Generalized Poincare inequality). For 

any \tp\ < tt there are constants c > and C < oo such 
that for all subsets O C K and all functions f on the 
torus K, the following estimate holds: 



1/2 



(25) 



Here |fi c | is the volume of fl c — K.\Q, the complement 
of in /C. 

Proof. We shall derive ( p5| ) from a special form of this 
inequality that holds for all functions that are orthogonal 
to the constant function. Namely, for any positive a < 
2/3 and some constants c > and C < oo (depending 
only on a and \tp\ < tt) we claim that 



T° ut = /i f dX f drxKVi+ip/LjttCn.X)! 2 , 



(21) 



2 i „ 



> 



L 2 



C 



m 

|/C| 



V' t llL 2 (K) ! 



(26) 



and 



i=\f dx/ dn^^Cln-rj-Dl^n.x) 



(22) 



provided (1, = 0. (Remark: Eq. ( |2^ ) holds also for 
a — 2/3, but -the proof is slightly more complicated in 
that case. SeeEJ.) If (p6| ) is known the derivation of j25| ) 
is easy: For any /, the function h = f — L~ 3 (1, f)ic is 
orthogonal to 1. Moreover, 



Here 

fi x = {n : |ri -r 3 -| < R for some j > 2} (23) 

is the complement of fix, and the diamagnetic inequality 
|(V + itp/L)f{r)\ 2 > |V|/(r)|| 2 has been used. The proof 
is completed by using the results ofEj ando (see alsoEl) 
which tell us that for e = Y 1 / 17 and R = aY- 5 / 17 

eT+(l- e)(T hl + I) > (l - const. Y 1/17 )4TTfipa (24) 

as long as N > Y~ x / 17 . (This estimate is highly non- 
trivial. Among several-other things it uses a generaliza- 
tion of Dyson's lemmaE3.) □ 



= l|V v /i]]| 2(x;) 



rJ L-I<1) _ T2 3/2 >/)k| 2 I 1 



|v^/|| 2 



+2|Re(i- 3 / 2 ,/) K (V v /,L- 3 / 2 }^ 



< l|V./" 2 



v>J \\L 2 (n) 



%\{L-* /2 ,f)K? 



-f(L\\V v f\\l HK) 



ii 



li 2 (/c) 



|/C| 



|O c | 

FT 



1/2 



(27) 
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and 



L 2 



\V v h\ 



< 



]2q) gives 



Setting a = J, using M » v'-n^ijy^ 
last term in (Eq) and combining ( |2q ) , (E7h and 

(||) with c = M + 

We now turn to the proof of (p6|). For simplicity we 
set L = 1. The general case follows by scaling. Assume 
that (pf ) is false. Then there exist sequences of con- 
stants C n — * oo, functions h n with ||/i n || i 2( K ) = 1 and 
(1, h n )tc — 0, and domains Q n C /C such that 



lim 



{l|V v Mi 2 (n„) + CnKHlV^M!^)} < ¥> 2 



(29) 



We shall show that this leads to a contradiction. 

Since the sequence h n is bounded in L 2 (tC) it has a sub- 
sequence, denoted again by h„, that converges weakly 
to some h <E L 2 (/C) (i.e., {g,h n )x. — > (g,h)ic for all 
.9 G L 2 (K,)). Moreover, by Holder's inequality the L P (Q^) 
norm || V^/i^Hz^c j = (/ nc M r )l p dr) 1/p is bounded by 

|^| Q/2 ||V y /i„|| i2(K) for p = 2/(a + !)• From @ we 
conclude that || V ¥ ,/i n ||ip(f2c ) is bounded and also that 
llV^/inllij.^) < \\V<phn\\L"{a n ) is bounded. Altogether, 
V v h n is bounded in L P (JC), and by passing to a further 
subsequence if necessary, we can therefore assume that 
V(o/i n converges weakly in L P (!C). The same applies to 
Vh n . Since p = 2/(a + 1) with a < r-?/3 the hypotheses 
of the Rellich-Kondrashov TheoremLD (Thm 8.9) are ful- 
filled and consequently h n converges strongly in L 2 (IC) to 
h (i.e., ||/i — h n \\ L 2^ — > 0). We shall now show that 



liminf ||V v /i n ||£ 2(n } > ||V y /i||^ 2(K) 



(30) 



This will complete the proof because the h n are normal- 
ized and orthogonal to 1 and the same holds for h by 
strong convergence. Hence the right side of |m]) is nec- 
essarily > (p 2 , since for \<p\ < n the lowest eigenvalue 
of — V 2 , with constant eigenfunction, is non-degenerate. 
This contradicts (p9|). 

Eq. ( |30"| ) is essentially a consequence of the weak lower 
semicontinuity of the L 2 norm, but the dependence 
on fl n leads to a slight complication. First, Eq. ( |29| ) 
and C n — ► oo clearly imply that |fi£j — > 0, because 
II ^v^«lli 2 (K) > </? 2 - By choosing a subsequence we 
may assume that ^ n < oo. For some fixed AT 
let f2jv = £ \ U„>jyf2°. Then Ojy C Q„ for n > N. 
Since || V v /i„|| 2 2( - \ is bounded, W v h n is also bounded 

in i 2 (f2jv) and a subsequence of it converges weakly in 
L 2 (D,m) to V v h. Hence 



liminf || V v h n \\ 2 L 2, nn) 



> liminf ||V<Ji n | 



IIV„/i 



(31) 



Since C f2jv+i and Un&n = (up to a set of 
measure zero), we can now let N — > oo on the right 
side of ( |3i"| ) . By monotone convergence this converges to 

j) which, as remarked above, 



I V(p/i||2,2(x;)- This proves 



— ||/-L- j (l,/) K ||i 2(K) . (28) contradicts @£ 



□ 



We now are able to finish the proof of Theorem [l]. From 
Lemmas I and | we infer that, for any symmetric \& with 
(ty, = 1 and for N large enough, 

±(*,H' N *)(l- const. Y 1 ' 17 )' 1 



> AiTfipa + /i—^ 



-^ 1/17 (^ + ^<^E i (v. 



-72 / dX 

Jk-n-i 



dri 



JK 

-L- 3 [J K dr^(r,X)] 



*(r ls X) 



(32) 



where we used that |ft c | < ^TVi? 3 = const. L 3 y 2 / 17 . 
From this we can infer two things. First, since the kinetic 
energy, divided by N, is certainly bounded independent 
of N, as the upper bound shows, we get that 



liminf — — — - — > Anapa 



/J 



L 2 



(33) 



for any \cp\ < ir. By continuity this holds also for \ip\ = w, 
proving (g). (To be precise, E /N — fitp 2 L~ 2 is concave 
in (p, and therefore stays concave, and in particular con- 
tinuous, in the limit N — > oo.) Secondly, since the upper 
and the lower bound to Eq agree in the limit considered, 
the positive last term in (^2|) has to vanish in the limit. 
I.e., we get that for the ground state wave function 
oiH' N 

lim / dX [ dri *o(ri,X) 

N^oo J k n-i J K 

-L- 3 [Xcdr* (r,X)]| 2 =0. (34) 
This proves (|^), since 

:1_ iV^/ ^ r > r ') drdr ' • ( 35 ) 

•J /C x K, 



dX / dri 

/C"- 1 JK 



and thereJbre N- 1 (L- 3 / 2 \~/ N \L- 3 / 2 ) -> 1. As ex- 
plained iioa this suffices for the convergence -/V _1 7jv — > 
L _3 / 2 )(i _3 / 2 | in trace class norm. □ 



IV. CONCLUSIONS 

We have shown that a Bose gas with short range, re- 
pulsive interactions is both a 100% superfluid and also 
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100% Bose-Einstein condensed in its ground state in the 
Gross-Pitaevskii limit where the parameter Na/L is kept 
fixed as N — > oo. This is a simultaneous large N and low 
density limit, because the dimcnsionless density param- 
eter pa 3 is here proportional to 1/N 2 . If pa 3 is not zero, 
but small, a depletion of the Bose-Einsteinj-condensate 
of the order (pa 3 ) 1 / 2 is expected (see, e.g.,E3). Never- 
theless, complete superfluidity in the ground state, e.g. 
of Helium 4, is experimentally observed. It is an inter- 
esting open problem to deduce this property rigorously 
from first principles. In the case of a one-dimensional 
hard-core Bose gas superfluidity in the ground state is 



easy to show, but nevertheless there is no Bose-Einstein 
condensation at all, not even in the ground stateaE3. 
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